On hadron production in Pb— Pb collisions at 158A GeV 
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A Fermi statistical model analysis of hadron abundances and spectra obtained in several relativistic 
heavy ion collision experiments is utilized to characterize a particle source. Properties consistent 
with a disintegrating, hadron evaporating, deconfined quark-gluon plasma phase fireball are ob- 
tained, with a baryochemical potential /is = 200-210 MeV, and a temperature T/ ~ 140-150 MeV, 
significantly below previous expectations. 



Discovery and study of quark-gluon plasma (QGP), a 
state consisting of mobile, color charged quarks and glu- 
ons, is the objective of the relativistic heavy ion research 
program underway at Brookhaven National Laboratory, 
New York and at CERN, Geneva ||l|. Thermalization of 
the constituents of the deconfined phase created in high 
energy large nuclei collisions is a well working hypothesis, 
as we shall see. The multi-particle production processes 
in 158yl GeV Pb-Pb coUisions carried out at CERN-SPS 
will be analyzed in this paper, using the principles of the 
statistical Fermi model |2j : strongly interacting particles 
are produced with a probability commensurate with the 
size of accessible phase space. Since the last comprehen- 
sive review of such analysis has appeared the Pb- 
beam experimental results became available, and model 
improvements have occurred: we implement refinements 
in the phase space weights that allow a full character- 
ization of the chemical non-equilibria with respect to 
strange and light quark flavor abundances |^J^. Con- 
sideration of the light quark chemical non-equilibrium 
is necessary in order to arrive at a consistent interpre- 
tation of the experimental results of both the wide ac- 
ceptance NA49-experiment ||^-[lo| and central rapidity 
(multi)strange (anti)baryon WA97-experiment |pl]-p^. 

We further consider here the influence of collective 
matter flow on to -particle spectra and particle multi- 
plicities obtained in a limited phase space domain. The 
different flow schemes have been described before [ p^ . 
We adopt a radial expansion model and consider the 
causally disconnected domains of the dense matter fire- 
ball to be synchronized by the instant of collision. We 
subsume that the particle (chemical) freeze-out occurs 
at the surface of the fireball, simultaneously in the CM 
frame, but not necessarily within a short instant of CM- 
time. Properties of the dense fireball as determined in 
this approach offer clear evidence that a QGP disinte- 
grates at Tf ~ 144 MeV, corresponding to energy den- 
sity e = 0(0.5) GeV/fm^ |l|]. Our initial chemical 
non-equilibrium results without How have been sugges- 
tive that this is the case |0 , showing a reduction of the 
chemical freeze-out temperature from Tf = 180 MeV ; 
an earlier analysis could not exclude yet higher hadron 



formation temperature of 270 MeV . 

The here developed model offers a natural understand- 
ing of the systematic behavior of the TO_L-slopes which 
differs from other interpretations. The near equality of 
(inverse) slopes of nearly all strange baryons and an- 
tibaryons arises here by means of the sudden hadroniza- 
tion at the surface of an exploding QGP fireball. In the 
hadron based microscopic simulations this behavior of 
TOj_-slopes can also arise allowing for particle-dependent 
freeze-out times 

In the analysis of hadron spectra we employ meth- 
ods developed in analysis of the lighter 200A GeV S- 
Au/W/Pb system where the description of the phase 
space accessible to a hadronic particle in terms of the pa- 
rameters we employ is given. Even though we use six pa- 
rameters to characterize the hadron phase space at chem- 
ical freeze-out, there are only two truly unknown prop- 
erties: the chemical freeze-out temperature Tf and light 
quark fugacity A, (or equivalently, the baryochemical po- 
tential /is = 3T/lnAg) - we recall that the parameters 
7i, i = q,s controls overall abundance of quark pairs, 
while Xi controls the difference between quarks and anti- 
quarks of given flavor. The four other parameters are 
not arbitrary, and we could have used their tacit and/or 
computed values: 

1) the strange quark fugacity As is usually fixed by the 
requirement that strangeness balances (s — s) = [||. 
The Coulomb distortion of the strange quark phase space 
plays an important role in the understanding of this con- 
straint for Pb-Pb coUisions, see Eq. (|^) p6[ ; 

2) strange quark phase space occupancy 7^ can be com- 
puted within the established kinetic theory framework 
for strangeness production 20 2l|]; 

3) the tacitly assumed equilibrium phase space occu- 
pancy of light quarks 7^ = 1 ; and 

4) assumed collective expansion to proceed at the rela- 
tivistic sound velocity, Vc = l/\/3 pl[ |. 

However, the rich particle data basis allows us to find 
from experiment the actual values of these four param- 
eters, allowing to confront the theoretical results and/or 
hypothesis with experiment. 

The value of As we obtain from the strangeness con- 
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servation condition (s - 
approximation llq]: 



: in QGP is, to a very good 



A,, = As A 
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Aq < 1 expresses the Coulomb deformation of strange 
quark phase space. This effect is relevant in central Pb- 
Pb interactions, but not in S-Au/W/Pb reactions. Aq is 
not a fugacity that can be adjusted to satisfy a chemical 
condition, since consideration of A^, i = u,d,s exhausts 
all available chemical balance conditions for the abun- 
dances of hadronic particles. The subscript i?/ in Eq. (|^) 
reminds us that the classically allowed region within the 
dense matter fireball is included in the integration over 
the level density. Choosing Rf ^ 8fni, T = 140 MeV, 
nis = 200 MeV (value of 7^ is practically irrelevant), for 
Zf = 150 the value is As = 1.10. 

In order to interpret particle abundances measured in 
a restricted phase space domain, we study abundance ra- 
tios involving what we call compatible hadrons: these 
are particles likely to be impacted in a similar fashion by 
the not well understood collective flow dynamics in the 
fireball. The available particle yields are listed in table |, 
top section from the experiment WA97 for pj^ > 0.7 GeV 
within a narrow Ay = 0.5 central rapidity window. Fur- 
ther below are shown results from the large acceptance 
experiment NA49, extrapolated to full 4tt phase space 
coverage. There are 15 experimental results. The to- 
tal error x| = " RLp)Vi^RLp? for the four 

theoretical columns is shown at the bottom of table ffl 
along with the number of data points W, parameters 
'p' used and (algebraic) redundancies 'r' connecting the 
experimental results. For r ^ it is more appropriate to 
quote the total x^, with a initial qualitative statistical 
relevance condition Xt/(-^^p) < 1- The four theoretical 
columns refer to results with collective velocity Vc (sub- 
script v) or without {vc — 0). We consider data including 
'All' data points, and also analyze data excluding from 
analysis four Q, particle ratios, see columns marked 
'No-f2. Only in letter case we obtain a highly relevant 
data description. Thus to describe the 17, yields we 
need an additional particle production mechanism be- 
yond the statistical Fermi model. We noted the special 
role of these particles, despite bad statistics, already in 
the analysis of the S-induced reactions 

Considering results obtained with and without flow 
reveals that the presence of the parameter Vc already 
when dealing only with particle abundances improves 
our ability to describe the data. However, m± spec- 
tra offer another independent measure of the collective 
flow Vc- for a given pair of values Tf and Vc we eval- 
uate the resulting particle spectrum and analyze it 
using the spectral shape and kinematic cuts employed 
by the experimental groups. To find the best values 
we consider just one 'mean' strange baryon experimen- 
tal value fl"" = 260 ± 10, since within the error the 
high m± strange (anti)baryon inverse slopes are over- 



lapping. Thus when considering Vc along with T± we 
have one parameter and one data point more. Once we 
find best values of Tf and Vc, we study again the in- 
verse slopes of individual particle spectra. We obtain 
an acceptable agreement with the experimental as 
shown in left section of table || . For comparison, we have 
also considered in the same framework the S-induced re- 
actions, and the right section of table |l| show a good 
agreement with the WA85 experimental data |25j. We 
used here as the 'mean' experimental slope data point 
= 235 ± 10. We can see that within a significantly 
smaller error bar, we obtained an accurate description 
of the m^-slope data. This analysis implies that the ki- 
netic freeze-out, where elastic particle-particle collisions 
cease, cannot be occurring at a condition very different 
from the chemical freeze-out. However, one pion HBT 
analysis at p± < 0.5 GeV suggests kinetic pion freeze- 
out at about T^ ~ 120 MeV |£6[. A possible explanation 
of why here considered p± > 0.7 GeV particles are not 
subject to a greater spectral deformation after chemical 
freeze-out, is that they escape before the bulk of softer 
hadronic particles is formed. 

The six statistical parameters describing the particle 
abundances are shown in the top section of table III, for 
both Pb- and S-induced reactions The errors shown 
are one standard deviation errors arising from the propa- 
gation of the experimental measurement error, but apply 
only when the theoretical model describes the data well, 
as is the case here, see the header of each column — 
note that for the S-induced reactions (see [||) the num- 
ber of redundancies is large since same data comprising 
different kinematic cuts is included in the analysis. We 
note the interesting result that within error the freeze- 
out temperature Tf seen in table III, is the same for 
both the S- and Pb-induced reactions. The collective 
velocity rises from — 0.5c to v^^ ~ = 0.58. 

We then show the light quark fugacity Xq, and note 
fi^^ = 203 ± 5 > /i| = 178 ± 5 MeV. As in S-induced 
reactions where As = 1, now in Pb-induced reactions, 
a value X^^ ~ 1.1 characteristic for a source of freely 
movable strange quarks with balancing strangeness, i.e., 
As = 1 is obtained, sec Eg. (|l|). 

7g > 1 seen in table |l| implies that there is phase 
space over- abundance of light quarks, to which, e.g., 
gluon fragmentation at QGP breakup prior to hadron 
formation contributes. 7^ assumes in our data analysis a 
value near to where pions could begin to condense [ p7| , 
7g = 7^ = e™'"/^^^ . We found studying the ratio / B 
separately from other experimental results that the value 
of 7q ~ 7^ is fixed consistently and independently both, 
by the negative hadron (ft.~), and the strange hadron 
yields. The unphysical range 7^ > 7^ can arise, since up 
to this point wc use only a first quantum (Bose/Fermi) 
correction. However, when Bose distribution for pions is 
implemented, which requires the constraint 7^ < 7^, we 
obtain practically the same results, as shown in second 
column of table III. Here we allowed only 4 free param- 
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and the strangeness conser- 
We then show in table [11 



eters, i.e. we set = 7^ 

vation constraint fixes As 

the ratio js/jq, which corresponds (approximately) to 
the parameter 7^ when 7^ = 1 had been assumed. We 
note that 7^'^ > 1. This strangeness over-saturation 
effect could arise from the effect of gluon fragmenta- 
tion combined with early chemical equilibration in QGP, 
7s (t < tf) ~ 1. The ensuing rapid expansion preserves 
this high strangeness yield, and thus we find the result 
7s > 1 , as is shown in figure 33 in |^ . 

We show in the bottom section of table III the 
energy and entropy content per baryon, and specific 
anti-strangeness content, along with specific strangeness 
asymmetry of the hadronic particles emitted. The energy 
per baryon seen in the emitted hadrons is nearly equal to 
the available specific energy of the collision (8.6 GeV for 
Pb-Pb, 8.8-9 GeV for S-Au/W/Pb). This imphes that 
the fraction of energy deposited in the central fireball 
must be nearly the same as the fraction of baryon num- 
ber. The small reduction of the specific entropy in Pb- 
Pb compared to the lighter S-Au/W/Pb system maybe 
driven by the greater baryon stopping in the larger sys- 
tem, also seen in the smaller energy per baryon content. 
Both collision systems freeze out at energy per unit of 
entropy E/S = 0.185 GeV. There is a loose relation of 
this universality in the chemical freeze-out condition with 
the suggestion made recently that particle freeze-out oc- 
curs at a fixed energy per baryon for all physical systems 
p8[ , since the entropy content is related to particle multi- 
plicity. The overall high specific entropy content we find 
agrees well with the entropy content evaluation made ear- 
lier Hi for the S-W case. 

Inspecting figure 38 in we see that the specific 
yield of strangeness we expect from the kinetic theory 
in QGP is at the level of 0.75 per baryon, in agreement 



with the results of present analysis shown in table 111 



This high strangeness yield leads to the enhancement of 
multi-strange (anti)baryons, which are viewed as impor- 
tant hadronic signals of QGP phenomena and a se- 
ries of recent experimental analysis has carefully demon- 
strated comparing p-A with A-A results that there is 
quite significant enhancement [|l3|,^, as has also been 
noted before by the experiment NA35 ]3^ ]. 

The strangeness imbalance seen in the asymmetrical S- 
Au/W/Pb system (bottom of table III) could be a real 
effect arising from hadron phase space properties. How- 
ever, this result also reminds us that though the statis- 
tical errors are very small, there could be a considerable 
systematic error due to presence of other contributing 
particle production mechanisms. Indeed, we do not offer 
here a consistent understanding of the f7, f2 yields which 
are higher than we can describe. We have explored addi- 
tional microscopic mechanisms. Since the missing fi, 
yields are proportional (13%) to the S, S yield, we have 
tested the hypothesis of string fragmentation contribu- 
tion in the final state, which introduces just the needed 
'shadow' of the S, S in the O, abundances. While this 
works for J7, f2, we find that this mechanism is not com- 



patible with the other particle abundances. We have also 
explored the possibility that unknown f2*, f2* resonances 
contribute to the fi, Q. yield, but this hypothesis is ruled 
out since the missing yield is clearly baryon-antibaryon 
asymmetric. Thus though we reached here a very good 
understanding of other hadronic particle yields and spec- 
tra, the rarely produced but greatly enhanced fi, U, must 
arise in a more complex hadronization pattern. 

We have presented a comprehensive analysis of hadron 
abundances and mj^-spectra observed in Pb-Pb 158^ 
GeV interactions within the statistical Fermi model with 
chemical non-equilibrium of strange and non-strange 
hadronic particles. The key results we obtained are: 
As = 1 for S and Pb colhsions ; 7^^ > 1, 7, > 1; 
S/B ~ 40; s/B ~ 0.75; all in a remarkable agreement 
with the properties of a deconfined QGP source hadroniz- 
ing without chemical re-equilibration, and expanding not 
faster than the sound velocity of quark matter. The uni- 
versality of the physical properties at chemical freeze-out 
for S- and Pb-induced reactions points to a common na- 
ture of the primordial source of hadronic particles in both 
systems. The difference in spectra between the two sys- 
tems arises in our analysis from the difference in the col- 
lective surface explosion velocity, which for larger system 
is higher, having more time to develop. Among other 
interesting results which also verify the consistency of 
our approach we recall: good balancing of strangeness 
(s — s) =0 in the Pb-Pb case; increase of the bary- 
ochemical potential as the collision system grows; en- 
ergy per baryon near to the value expected if energy and 
baryon number deposition in the fireball are similar. We 
note that given the magnitude of 7^ and low chemical 
freeze-out temperature, most (75%) of all final state pi- 
ous are directly produced, and not resonance decay prod- 
ucts. Our results differ significantly from an earlier anal- 
ysis regarding the temperature at which hadron forma- 
tion occurs. Reduction to Tf — 140-145 MeV becomes 
possible since we allow for the chemical non-equilibrium 
and collective flow, and only with these improvements 
in analysis our description acquires convincing statistical 
significance, which e.g. a hadronic gas scenario does not 
offer Q. Because we consider flow effects, we can ad- 
dress the central rapidity data of the WA97 experiment 
at the required level of precision, showing the consistency 
in these results with the NA49 data discussed earlier . 

In our opinion, the only consistent interpretation of 
the experimental results analyzed here is that hadronic 
particles seen at 158^ GeV nuclear collisions at CERN- 
SPS are formed directly in hadronization of an exploding 
deconflned phase of hadronic matter, and that these par- 
ticles do not undergo a chemical re-equilibration after 
they have been produced. 
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TABLE I. WA97 (top) and NA49 (bottom) Pb-Pb 158A 
GeV particle ratios and our theoretical results, see text for 
explanation. 



Ratios Ref. Exp. Data 



H/A 
E/A 
A/A 

n/E 

Ti/E 
Ti/Q. 
(n+n) 

(5+S) 



(S+5) 
(A+A) 

K+/K' 

Pjp 
A/p 
K°/B 



|12| 
|12| 
12 

in 

|12 

111 
12 

fl 



0.099 ± 0.008 
0.203 ± 0.024 
0.124 ± 0.013 
0.255 ± 0.025 
0.192 ± 0.024 
0.27 ± 0.06 
0.38 ± 0.10 

0.20 ± 0.03 



22 

1 

1^1 
ll 
|24 
|23 
23 



0.13 ± 0.03 
11.9 ± 1.5 
1.80± 0.10 
18.1 ±4. 
3. ± 1. 
0.183 ± 0.027 
1.83 ± 0.2 



N;p;r 



AU 



AU|„ 



0.107 

0.216 

0.121 

0.246 

0.192 

0.40 

0.51 

0.23 



0.110 

0.195 

0.128 

0.225 

0.190 

0.40 

0.47 

0.23 



0.109 

16.1 

1.62 

16.7 

0.65 

0.242 

1.27 



0.111 

15.1 

1.56 

15.3 

1.29 

0.281 

1.55 



19 

15;5;4 



18 
16;( 



No-Q. No-QL 



0.095 

0.206 

0.120 

0.260 

0.078* 

0.17* 

0.57* 

0.10* 



0.102 

0.210 

0.123 

0.252 

0.077* 

0.18* 

0.60* 

0.10* 



0.107 

9.89 

1.76 

17.3 

2.68 

0.194 

1.80 



0.114 

12.9 

1.87 

17.4 

2.02 

0.201 

1.83 



2.1 

11;5;2 



1.8 

12;6;2 



TABLE H. Experimental and theoretical m± spectra in- 
verse slopes Tth. Left Pb-Pb results from experiment NA49 
[ p^ for kaons and from experiment WA97 for baryons; 
right S-W resuhs from WA85 S. 





Tl^ [MeV] 


Tth" [MeV] 


Tl [McV] 


[MeV] 




223 ± 13 


241 


219 ± 5 


215 


rpA 


291 ± 18 


280 


233 ± 3 


236 


rj-iA 


280 ± 20 


280 


232 ± 7 


236 


T- 


289 ± 12 


298 


244 ± 12 


246 




269 ± 22 


298 


238 ± 16 


246 



TABLE III. Top section; statistical parameters, and 
their Xtj which best describe the experimental results for 
Pb-Pb data, and in last column for S-Au/W/Pb data pre- 
sented in Ref. . Bottom section: specific energy, entropy, 
anti-strangeness, net strangeness of the full hadron phase 
space characterized by these statistical parameters. In the 
middle column we fix As by requirement of strangeness con- 
servation and choose 7q = 7^, the pion condensation point. 





Pb-No-n|„ 


Pb-No-n|: 


S-No-nj„ 


Xt; N;p;r 


1.8; 12; 6; 2 


4.2; 12; 4; 2 


6.2; 16; 6; 6 


Tl [MeV] 


144 ± 2 


145 ± 2 


144 ± 2 


Vc 


0.58 ± 0.04 


0.54 ± 0.025 


0.49 ± 0.02 


\ 


1.60 ± 0.02 


1.605 ± 0.025 


1.51 ± 0.02 


\s 


1.10 ± 0.02 


1.11* 


1.00 ± 0.02 


Iq 


1.7 ± 0.5 




1.41 ± 0.08 


Ishq 


0.86 ± 0.05 


0.78 ± 0.05 


0.69 ± 0.03 


Ej/B 


7.0 ± 0.5 


7.4 ± 0.5 


8.2 ± 0.5 


Sf/B 


38 ± 3 


40 ± 3 


44 ± 3 


Sf/B 


0.78 ± 0.04 


0.70 ± 0.05 


0.73 ± 0.05 


(sf - Sf)/B 


0.01 ± 0.01 


0* 


0.17 ± 0.02 
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